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Extending the Eigenvector 1 Space to the Optical Variability of 

Quasars 

Yufeng Mao 1,2 , Jing Wang 1 , & Jianyan Wei 1 
ABSTRACT 

We introduce a new physical parameter, the optical variability amplitude, to 
the well-established Eigenvector 1 space of quasars and test a sample of long-term 
B-band light curves of 42 PG quasars monitored by Giveon et al. (1999). We 
find that the optical variability amplitude strongly correlates with the intensity 
ratio of Fe II to H/?, H/3 width and peak luminosity at 5007A. We briefly discuss 
the physical meaning of our findings and suggest that the Eddington ratio may 
be a key factor in determining a quasar's variability. 

Subject headings: galaxies: active — quasars: general 



1. Introduction 

During the last few decades, quasars have been monitored in multi-wavelength observa- 
tions, from radio to X-rays, by many research programs. Our understanding of the central 
engine of quasars has been greatly improved by the correlations found between the variabil- 
ity properties and other observational parameters (e.g., luminosity, redshift, rest-wavelength, 
timescales and emission-line width). An anticorrelation between the amplitude of variability 
and luminosity was reported by Pica & Smith (1983), and confirmed by many subsequent 
studies (e.g., Cristiani et al. 1990; Hook et al. 1994; Cristiani et al. 1997; Giveon et al. 
1999 (hereafter G99); Garcia et al. 1999; Vanden Berk et al. 2004), although there have 
also been reports to the contrary (e.g., Trevese et al. 1989; Giallongo et al. 1991; Cimatti 
et al. 1993). The relationship between the variability amplitude and redshift was discussed 
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in many studies. Some authors found that the variability amplitude is anticorrelated with 
redshift (e.g., Barbieri et al. 1983; Cristiani et al. 1990; Hook et al. 1994; Cristiani et al. 
1996), while others reported an opposite trend (e.g., Giallongo et al. 1991; Treverse et al. 
1994; Cid Fernandes et al. 1996). Moreover, it was found that the optical spectra usually 
become harder as quasars turn brighter (Cutri et al. 1985; G99). Kollatschny et al. (2006) 
examined the variability properties of a sample of 10 Palomar-Green quasars with the line 
width of H/5 larger than 5000 km s -1 . They found a marginal correlation between the optical 
continuum variability amplitude and H/5 line width, which provides useful information for 
understanding the structure of BLR. 

On the other hand, Boroson & Green (1992; hereafter BG92) examined a sample of 87 
bright low redshift PG quasars and found that most of the variance is connected to two sets of 
correlations, which were then defined as Eigenvector 1 (hereafter El) space and Eigenvector 
2 space. The El space is dominated by the strong anticorrelation between the optical Fell 
and [OIII], and Eigenvector 2 by the correlation between the optical luminosity and Hell 
equivalent width (see Sulentic et al. 2000a for a review). By calculating the virial black hole 
mass (Mbh) using the well established empirical R — L relationship (e.g., Kaspi et al. 2000), 
Boroson (2002) suggested that the El space is dominantly driven by the Eddington ratio 
L/LEdd, and the Eigenvector 2 by Mbh (see also in Sulentic et al. 2000b). 

In this paper, we will investigate whether the quasar's optical variability amplitude is 
related with the El space. We test a sample of 42 PG quasars. The paper is structured 
as follows. The sample selection and our analysis are described in Section 2. The results 
are presented in Section 3, and discussions in Section 4. Throughout this paper, the ACDM 
cosmology (Spergel et al. 2003) with VL m = 0.3, = 0.7, and H = 70 km s _1 Mpc -1 are 
adopted. 

2. Sample and Analysis 

We searched the literature for suitable results of quasars' optical variability. In partic- 
ular, the light curves should be well sampled in the temporal domain, with adequate total 
observation time and sampling interval. It is found that the G99 subset of the PG quasars 
is well suited for our purpose, because the optically selected PG quasars are not only nearly 
statistically complete but also studied comprehensively in multi-wavelength observations. 
G99 monitored 42 nearby (z < 0.4), bright (B < 16mag) PG quasars in the B and R bands 
for a seven-year period at Wise Observatory. The typical temporal sampling interval was 
40 days, and the objects were observed at 30-60 epochs with a photometry uncertainty of 
~ O.Olmag. All the objects showed intrinsic rms variability amplitudes of 5% < cr B <34% 
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and 4% < a R <26% . 

Our sample is listed in Table I. Figure la shows the distribution of the redshifts for the 
sample listed in G99, see also Column (2) of Table 1. The redshifts of these PG quasars are 
mainly less than 0.2. Column (5) lists the variability amplitude in magnitude for each object. 
Furthermore, the variability amplitude can be statistically estimated in various ways. We 
refer the readers to G99 for a brief comment on these different methods. In the current study, 
the variability is defined as the median value of all possible magnitude difference of a light 
curve, simply because the median value is a relatively robust estimation, i.e., not strongly 
affected by the outliers (e.g., Hook et al. 1994, Netzer et al. 1996). Only the variability in 
the B band is considered in the subsequent analysis since the variability is more significant 
in the B band than in the R band. G99 examined the correlation between the variability 
amplitude defined by magnitude and quasar spectral properties defining the El. However, 
no significant correlations were found by them. It should be pointed out that the variability 
amplitude defined in magnitude represents a relative change in luminosity. For a constant 
change defined in luminosity, a small (large) change in magnitude could be simply caused 
if the object is (less) luminous. In order to overcome this problem, the median of absolute 
change in luminosity is used in this paper, that is 

AL = L bol x (1 - i -°- 4 x|"^"( AB )l) (!) 

where median(AB) is the median of variability in the B band and Lboi the fiducial luminosity 
of the quasar. Note that AL is always positive according to this definition. AL represents 
the characteristic variability of each quasar in the absolute luminosity change, which directly 
reflects the absolute change of the amount of the fueling gas. Regarding two quasars with 
the same change in the median(AB), the one with a larger bolometric luminosity would 
have a larger AL. 

The bolometric luminosity can be estimated in two ways. First, we use an widely 
accepted empirical relationship 

L bo i = 9AL A (5100l) (2) 

as given by Kaspi et al. (2000), where the luminosities at the rest-frame wavelength 5100A 
are adopted from Vestergaard & Peterson (2006). In another way, we also calculate the L bo i 
from the median apparent magnitudes in the B band given by G99, using the formula 

log(L bol /v B L VB ) = 0.80 - 0.067£ + 0.017£ 2 - 0.0023£ 3 (3) 

of Marconi et al. (2004), where C = (log L bo i -12). A comparison of between the values of 
Lboi obtained from Equation (2) and those from Equation (3) is made (see Figure 2). There 
is a strong correlation (with slope ~ 1) between the two sets of L bo i, indicating that the two 
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independent measurements are highly consistent with each other. So we adopt the value 
of Lboi obtained from Equation (2), and then calculate the AL using Equation (1). The 
calculated AL in the logarithm is shown in Column (6) of Table 1 for each quasar. 

3. Results 

The main goal of the present paper is to investigate whether the variability amplitude of 
quasars is related to El space. The El space is dominated by significant correlations between 
RFe (=FeII/H/3), FWHM(H/3), and [OIII] strength, and has been discussed by many authors 
(e.g. BG92; Xu et al. 2003; Grupe 2004; Sulentic et al. 2000a). We list the El parameters 
of our PG quasar sample in Table 1. Column (7) gives the FWHM of the broad component 
of E(3, Column (8) the ratio of the peak hight of [OIII] A5007 to that of H/3 (Peak A5007), 
Column (9) the ratio of the flux of Fell integrated in the rest frame wavelength range from 
A4434 to A4684 to that of R/3 (RFe), Column (10) the logarithm of R, i.e., the ratio of radio 
flux at 6 cm to optical flux density. All of these parameters are adopted from BG92. 

We then investigate the correlations between AL and El parameters in our sample. 
Figure 3 displays the correlations of AL versus FWHM(H/3), peak A5007, RFe, and the radio 
loudness log R, respectively. The spearman rank-ordered correlation coefficients r s of the four 
correlations are listed in Table 2, where P s is the probability of null correlation. Figure 3a 
shows a significant correlation between the AL and FWHM(H/9). A spearman rank-ordered 
analysis yields a correlation coefficient r s =0.450 with a significance level P s =0.004. This 
means that the quasars with larger widths of H(3 would have larger changes in luminosity. We 
also find a correlation between AL and peak A5007 (Fig. 3b, r s =0.445, P s =0.004). The anti- 
correlation between AL and RFe is plotted in Fig.3c. The calculated correlation coefficient 
is r s =-0.441, and the significance level P s =0.005. Kollatschny et al. (2006) did not find a 
significant correlation between the continuum variability amplitude at 5100A and radio power 
at 5GHz in their sample. In current studies, a correlation (r s =0.476, P s =0.0023) is identified 
between the radio loudness logP and AL, and is shown in Fig. 3d. Since logP = 1 is widely 
used to separate the radio-loud and radio-quiet quasars (Kellermann et al. 1989), the diagram 
shows that all the radio-loud quasars (logP > 1) have large optical variability amplitudes 
(log(AL) > 45), although the radio-quiet quasars (logP < 1) are nearly evenly distributed 
in terms of variability amplitude. The fact that radio-loud quasars have large variability 
amplitudes implies that the optical continuum of radio-loud quasars is contaminated by the 
high energy tail of their radio emissions, which boosts the variability amplitude because of 
the beam effect of the jets. 

The correlations found above suggest that the El space could be extended to include 
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the variability amplitude. This hypothesis can be verified by a principal component analysis 
(PCA) of our sample. The PCA is performed using the following 11 parameters, which are: 
log.R, the equivalent width of H/3, R(A5007), R(A4686), RFe, Peak(A5007), FWHM(H/3), H/3 
shift, B.(3 shape, H[3 asymmetry, and log(AL), each potentially providing unique information. 
Except for log(AL), the former ten parameters are directly collected from BG92. We refer 
the reader to BG92 for the definitions of these parameters. The PCA results are presented 
in Table 3, which lists the first four most significant eigenvectors. The second row shows the 
cumulative percentage of the variance. One can see that the first four eigenvectors together 
account for more than 70 percent of the variance, and that the first principal component 
dominates the observed properties of quasars. Similar to BG92, the El is dominated by 
the anticorrelation between the strength of Fell and Peak (A5007). In addition, our El is 
strongly effected by log(AL). It is clear that log(AL) has a projection of 0.54 on El, and 
0.57 on the Eigenvector 2. Although at first sight, the projection on Eigenvector 2 is larger 
than on El, taking the larger cumulative percentage of the El into account, we conclude 
that the El space can be extended to log(AL). 



4. Discussion 

4.1. Variability vs. M B h and Eddington ratio 

After extending the El space to the variability in amplitude, the dominant physical 
parameters are discussed in this section. Both of the black hole mass (M B h) and Eddington 
ratio (L/L Edd ) are believed to be the main parameters governing the observed properties 
in quasars. However, the luminosity variability defined by Equation (1) cannot be used to 
correlate with Mbh and L/LEdd because both Mbh and L/LEdd are estimated in terms of the 
continuum luminosity which is used to define the variability in luminosity. The variability 
in magnitude median(AB) is therefore used instead in the subsequent analysis. All Mbh 
values are adopted from Vestergaard & Peterson (2006). The distribution of Mbh is shown 
in Figure lb. One can see that the majority of log(MBH/M Q ) lie between 7.5 and 9.5. Lboi is 
then estimated using the Equation (2). The values of M B h and L/L E dd are listed in Column 
(3) and (4) in Table 1, respectively. 

The relation between the variability amplitude and Mbh or L/LEdd has been discussed 
recently. Contradictory results were, however, obtained by different authors. Wold et al. 
(2007) examined the relation between the quasar variability and black hole mass by studying 
the optical variability of ~100 quasars monitered by QUEST1 survey (Rengstorf et al. 2004). 
However, a correlation between the R-band variability and M B h was marginally obtained only 
when Mbh was averaged in several bins. Furthermore, they did not detect such a relation 
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in their PG quasar sample. In contrast, Wilhite et al. (2008) found a correlation between 
the variability and M B h in a sample of ~2500 quasars selected from the SDSS. In addition, 
they reproduced the well-known anti-correlation between the variability and luminosity. By 
combining the two relations, they suspected that L/L Edd is a possible driver for quasar 
variability in an indirect way. 

As Mbh (FWHM) 2 , a possible way to test the correlation between variability and 
Mbh is to search for the correlation between the variability and line width, both of which 
are independent observational parameters. The relation between the variability amplitude 
and line width was discussed in several papers. G99 found a marginal correlation between 
the variability amplitude (defined in magnitude) and width of the Hf3 emission line. Their 
possible, but unlikely explanation, for this trend is the contributions of the emission lines to 
the broad-band emission. Kollatschny et al. (2006) confirmed the results of G99 in a sample 
of 43 galaxies. In the above analysis, we find a significant correlation between log(AL) and 
the FWHM of H/3 (Fig. 3a). 

The median (AB) is plotted against M B h in Figure 4a. No significant correlation is, 
however, found between these two parameters (r s = 0.124, P s = 0.4267). Wold et al. (2007) 
did not detect a correlation between variability and Mbh in their PG quasar sub-sample, 
and our result confirms their conclusion. However, the current results make it difficult to 
understand why the correlation between the magnitude variability and the Mbh is not as 
good as expected. 

It is now generally believed that the El spaces is likely driven by L/L E dd (e.g., Boroson 
2002). The relation between the magnitude variability and L/LEdd is directly examined 
in Figure 4b. We find a significant anti-correlation between the median(AB) and L/LEdd 
(r s = -0.368, P s = 0.0012). Although our result agrees with Wilhite et al. (2008), caution 
must be made when explaining the median(AB)-L/LEdd correlation. Taking two quasars 
with the same Mbh, the more luminous one would have smaller variability in magnitude for 
a given change in luminosity (accretion rate). Meanwhile, the more luminous quasar would 
have a larger L/L E dd- That means that the median(AB)-L/L E dd correlation might be caused 
by an intrinsic relation in mathematics (i.e., the definition of magnitude) rather than in 
physics. A sample of light curves defined in flux or luminosity is therefore required to test 
the underlying physics of the relation. 
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4.2. Implications on variability mechanisms 

Several theoretical models have been proposed as mechanisms for variability in quasars, 
such as disk-instability (Kawaguchi et al. 1998), gravitational microlensing (Hawkings 1993, 
1996), and starburst (Terlevich et al. 1992). However, all of them are still far from clear. 

Hawkins (1993, 1996) explained the observed AGN variability by invoking gravitational 
microlensing. In this model, a quasar's light is lensed by a large population of compact 
bodies with planetary-mass. The microlensing model has two parameters: the Einstein 
radius of the lenses and their mean transverse velocity (Hawkins 2002). However, the two 
parameters generally can not be obtained observationally. In addition, microlensing events 
should be extremely rare at low redshift (Vanden Berk et al. 2004). This explanation can 
be easily excluded for the PG quasars studied in this paper. We find the quasar's variability 
is related with the El space, which strongly indicates that the variability must be caused by 
an intrinsic mechanism rather than an external one. 

Someone holds the idea that AGN variability might be caused by a series of discrete 
outbursts, such as supernova explosions (Aretxaga et al. 1997). However, this model can 
not explain the relationship between the luminosity and variability amplitude as argued by 
Pica & Smith (1983). Alternatively, Terlevich et al. (1992) explained the AGN variability 
as originating from the supernova remnants (SNRs) occurring in the innermost regions of 
AGNs. The long-term and short-term variability observed in AGNs could be explained by 
the long-term decay of the SNRs and cooling instability after the onset of the radiative phase, 
respectively. The cooling time before the radiative phase is ~0.6 year, and beyond this, the 
phase is reduced to ~6 days. 

The disk-instability model is much more popular than the other models. This model 
interprets the variability as an occasional flare event or blob formation caused by the insta- 
bility in the accretion disk. Kawaguchi et al. (1998) compared the logarithmic slopes of the 
structure function between the disk-instability model and star-burst model, and their obser- 
vation of quasar 0957+561 supports the disk-instability model. Vanden Berk et al. (2004) 
studied photometric variability of 25,000 quasars from SDSS, and found that their results 
favor the disk-instability model. However, it is still unclear how changes in the accretion rate 
or the resulting luminosity changes would propagate through the accretion disk. Recently, Li 
& Cao (2008) proposed a disk model, and suggested that the disk temperature change would 
lead to systematic spectral shape difference, which could explain the correlation between the 
variability and L/L E dd or M B h discovered by Wold et al. (2007) and Wilhite et al. (2008). 

Although our analysis can not discriminate which model, starburst or disk-instability, is 
favored for a quasar's variability, the extension of the El space to the variability implies that 
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either: 1) the stability of the accretion rate (or gas supply) changes along the El sequence 
in the disk-instability model, or 2) the intensity of star formation activity changes along the 
El sequence. 

5. Conclusions 

By studying the variability of the 42 PG quasars monitored by Giveon et al. (1999), 
using both direct correlations and PCA analysis, we find that the El space can be extended 
to include quasar's optical variability. The link between this variability and Eddington 
ratio/black hole mass is discussed, and we propose that the Eddington ratio may be a key 
factor in determining the variability of quasars. 
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of China. 
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Fig. 1. — Distribution of redshift (Panel A) and black hole mass (Panel B) for the 42 PG quasars 
studied in this paper. 
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Fig. 2. — Comparison between the value of Lboi calculated from Equation 2 and 3. The horizon 
axis represents the Lboi from Equation 2 (Kaspi et al. 2000), and the vertical axis represents the 
L bol from Equation 3 (Marconi et al. 2004). It is clear that the two independent measurements are 
consistent with each other. 
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(r s =0.445, P s =0.0044). (c) RFell (r s =-0.441, P s =0.0048). (d) log R(r s =0.476, P s =0.0023). 
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Table 2: Spearman Rank-Order Correlation Coefficient of the correlations shown in Figure 3 



A L vs 


FWHM(H/3) 


Peak A5007 


RFe 


Log R 


r s 


0.450 


0.445 


-0.441 


0.476 


Ps 


0.004 


0.004 


0.005 


0.002 



Table 3: Correlations of Eigenvectors with line and continuum properties 



Property 


Eigenvector 1 


Eigenvector 2 


Eigenvector 3 


Eigenvector 4 


Eigenvalue 


3.89 


1.87 


1.24 


1.03 


Cumulative 


35.3% 


52.3% 


63.6% 


73.0% 


logR 


0.68 


0.54 


-0.15 


0.03 


EW(H/3 b ) 


-0.05 


-0.60 


-0.55 


-0.08 


R (A5007) 


0.82 


0.04 


0.43 


-0.13 


R (A4686) 


-0.14 


-0.55 


0.44 


-0.23 


RFe 


-0.79 


0.45 


-0.03 


-0.04 


Peak (A5007) 


0.92 


-0.07 


0.24 


-0.12 


FWHM(H/3 b ) 


0.68 


-0.21 


-0.32 


0.34 


H/3 shift 


-0.10 


0.39 


0.46 


0.05 


H/3 shape 


0.07 


0.27 


-0.29 


-0.89 


H/3 asymm 


-0.70 


-0.33 


0.02 


0.15 


log (AL) 


0.54 


0.57 


-0.30 


0.10 



